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We report improved measurements of B to pseudoscalar-vector decays containing an ui meson 
in the final state. Our results are obtained from a data sample that contains 388 x 10 6 BB pairs 
accumulated at the T(4S) resonance, with the Belle detector at the KEKB asymmetric-energy 
e + e~ collider. We measure the following branching fractions: B(B + — > ujK + ) = [8.1 ± 0.6(stat.) ± 
0.6(syst.)] x 1(T 6 , B{B + -> luk + ) = [6.9 ± 0.6(stat.) ± 0.5(syst.)] x 1(T 6 , and B(B° -> luK°) = 

[4.4±°;®(stat.)±0.4(syst.)] X 10~ 6 . The partial width ratio ^^ZZk") = l-7±0.3(stat.) ±0.1(sys.). 
We also set the 90% confidence level upper limit B(B° — > loty ) < 2.0 x 10 -6 . In addition, we 
obtain the partial rate asymmetries Acp = 0.05jlo]o|(stat.) ± 0.01(syst.) for B + — > ujK + , and 
Acp = -0.02 ± 0.09(stat.) ± 0.01(syst.) for B + -* u-k + . 

PACS numbers: 11.30.Er,13.25.Hw 



Charmless hadronic B decays provide a rich ground to 
understand the dynamics of B meson decays and the ori- 
gin of CP violation. Two-body B decays with a vector 
meson and a pseudoscalar particle h (h is either a kaon 
or a pion) proceed through combinations of color-allowed 
(T u ) and color-suppressed (C u ) b — > u tree diagrams and 
b — ► s (P s ) or b — > d (Pd) penguin diagrams. For ex- 
ample, B + — > loK + proceeds through T u , C u and P s , 
B° — > uiK° proceeds through P s , B + — > ojit + proceeds 
through T Ui C u and Pj, and B° — > ujtt proceeds through 
C u and Pd- The QCD factorization (QCDF) approach 
suggests that the branching fractions of B + — ► ojh + 
decays are in the range 10~ 6 -10 -5 The predicted 
B(B° -> ujK°) and B{B+ -> loK+) are enhanced after 
including the rescattering effect from final-state inter- 
actions [2|]. The perturbative QCD (PQCD) approach, 
on the other hand, suggests that B(B + — > ujK + ) = 
3.22 x 1CT 6 , and B(B° uK a ) = 2.07 x 10~ 6 , and 
that B(B+ ojK+) = (10.6+5° 8 Ji+ 7 4 J ) ) x 10" 6 , and 

B(B° -» ojK°) = (9.8+®; ^9(^4.3)) x 10 ~ 6 after consider- 
ing ncxt-to-leading-order accuracy Due to the lack 
of T u and P s diagrams, B° — > ujtt is expected to be 
small 0, Q . Experimentally, clear signals have been ob- 
served for B+ -> u>K+, B+ -> W7T+, and B° -» uK° with 
similar branching fractions [1, @, 0, S| ■ However, experi- 
mental measurements to date are not yet precise enough 
for a quantitative confirmation of the pattern predicted 
by QCDF or PQCD. 



In this paper, we report improved measurements of 
branching fractions and partial rate asymmetries for B — ► 
tuh decays, where h can be a kaon or pion. The partial 
rate asymmetry (Acp) is measured for the charged B 
decays and defined to be 

_ T(B~ ^uh-)-r(B+^Ljh+) 
CP ~ T(B- -► uh~) + T(B+ wh+) ■ 1 ' 

These measurements are based on a data sample of 
388 x 10 6 BB pairs collected with the Belle detector at the 
KEKB Q asymmetric-energy e + e~ collider. They im- 
prove upon our previously published results [f| by a five- 
fold increase in statistics and supersede them. The Belle 
detector is a large-solid-angle magnetic spectrometer that 
consists of a silicon vertex detector (SVD), a 50-layer 
central drift chamber (CDC), an array of aerogel thresh- 
old Cerenkov counters (ACC), a barrel- like arrangement 
of timc-of- flight scintillation counters (TOF), and an 
electromagnetic calorimeter (ECL) comprised of CsI(Tl) 
crystals located inside a superconducting solenoid coil 
that provides a 1.5 T magnetic field. An iron flux-return 
located outside of the coil is instrumented to detect K\ 
mesons and to identify muons (KLM). The detector is 
described in detail elsewhere [10(. In August 2003, the 
three-layer SVD was replaced by a four-layer radiation 
tolerant device. The data sample for this analysis con- 
sists of 140 fb _1 of data recorded with a three-layer SVD 
(Set I) [3 and 217 fb" 1 recorded with a four-layer SVD 
(Set II) TU. 



Hadronic events are selected using criteria based on the 
charged track multiplicity and total visible energy, with 
an efficiency greater than 99% for generic BB events. 
All primary charged tracks must satisfy quality require- 
ments based on their impact parameters relative to the 
run-dependent interaction point (IP). For tracks from the 
candidate B mesons, their deviations from the IP posi- 
tion are required to be within ±0.1 cm in the transverse 
direction and ±3.0 cm in the longitudinal direction. 

Charged particle identification is performed using a K- 
tt likelihood ratio, TZk = C-k / (Ar + where Ck (At) 
is the likelihood for a charged particle to be a kaon (pion) 
based on information from the ACC, TOF and CDC. 
Charged tracks with TZk > 0.6 are identified as kaons, 
and tracks with TZk < 0.4 are identified as pions. With 
these criteria, kaons produced in B + — > ljK + decays are 
selected with an efficiency of 85%, while the correspond- 
ing rate of kaons that are misidentified as pions is 11%. 



On the other hand, pions from B + 



are selected 



with an efficiency of 90% and have a corresponding kaon 
misidentification rate of 9%. 

Candidate ir° mesons are reconstructed from pairs 
of photons with an invariant mass in the range of 
118 MeV/c 2 < Mp, < 150 MeV/c 2 (within ±2.5ct around 
the PDG value 12|). The cosine of the decay angle 
should satisfy |cos# 7 | < 0.9, where # 7 is the angle be- 
tween the photon decay axis and the negative of the 
laboratory frame direction in the 7T° rest frame. Can- 
didate Kg mesons are reconstructed using pairs of op- 
positely charged tracks that have an invariant mass in 
the range of 482 MeV/c 2 < M w+7r - < 514 McV/c 2 
(within ±4<7 around the PDG value [13 )■ The vertex 
of the K s candidate is required to be well reconstructed 
and displaced from the IP, and the K s momentum di- 
rection must be consistent with the K s flight direction. 
Candidate uj — > 7r + 7r~7r° decays are reconstructed from 
charged tracks with TZk < 0.9 (this requirement has a 
96% efficiency per track), and from 7r° candidates with 
e + e~ centcr-of-mass frame (CM) momenta greater than 
0.35 GeV/c 2 . Candidate lu mesons are required to have 
invariant masses within ±30 McV/c 2 (~ ±3er) of the 
PDG value 0. 

B meson candidates are formed by combining an to 
meson with either a kaon (K + , Kg) or a pion (tt + , 
7T°). Two kinematic variables are used to select B 
candidates: the beam-energy constrained mass M^ c = 

(£™ m ) 2 - |P£ M | 2 , and the energy difference AE = 

E% M -EP™, where Ep™ is the beam energy in the CM 
frame, and Pg M , Eg M are the momentum and energy, 
respectively, of the B candidate in the CM frame. Can- 
didates with Af bc > 5.2 GeV/c 2 and \AE\ < 0.25 GeV 
(\AE\ < 0.30 GeV for B° -> lutt ) are selected for further 
analysis. 

The dominant background arises from quark-antiquark 
continuum events (e + e~ — > qq, q = u,d,s,c). The contin- 



uum background is characterized by a jet-like structure, 
while the BB events have a more spherical distribution. 
The following event-shape variables calculated in the CM 
frame are employed to suppress the continuum. The 
thrust angle &t is defined as the angle between the thrust 
axis [l3| of the B candidate daughter particles and that 
of the rest of the particles in an event. Signal events are 
uniformly distributed in cos 8t , while continuum events 
are sharply peaked near cos 9t = ±1.0. Events with 
| cos 0t I < 0.9 are selected. 

A Fisher discriminant is formed by combining a set 
of modified Fox- Wolfram moments [lj] with the variable 
S± [l5| . This variable is the scalar sum of the transverse 
momenta of particles outside a 45° cone around the B 
thrust axis, divided by the scalar sum of their momenta. 
Here, the transverse momenta are calculated with respect 
to the thrust axis, and we do not include daughters of the 
B candidate. Further variables that have been found to 
separate signal from continuum background include: the 
cosine of the angle between the flight direction of the 
B candidate and the beam direction (cos6*b); the dis- 
tance along the beam direction between the B vertex 
and the vertex of the remaining particles in the event 
(Az)\ and the cosine of the helicity angle, defined as the 
angle between the direction opposite to the B flight di- 
rection in the lo rest frame and the direction normal to 
the plane defined by the three daughter pions of the to. 
The probability density functions (PDFs) for these three 
variables and the Fisher discriminant are obtained us- 
ing Monte Carlo (MC) simulation for signal events and 
sideband data (M\, c < 5.26 GeV/c 2 ) for qq backgrounds. 
These variables are combined to form a likelihood ratio 
^ = Cs+ S c - ' wnere £-s(qq) is the product of signal (qq) 
PDFs. 

Additional background discrimination is provided by 
the quality of the B flavor tagging of the accompany- 
ing B meson. We use the standard Belle B tagging 
package [3], which gives two outputs: a discrete vari- 
able (q) indicating the B flavor, and a dilution factor 
(r) ranging from zero for no flavor information to unity 
for unambiguous flavor assignment. We divide the data 
into six r regions. Continuum suppression is achieved 
by applying a mode-dependent requirement on TZ for 
events in each r region based on the figure-of- merit: 
•Ms/ y/N s + Nqq, where Af s is the expected number of sig- 
nal events estimated from simulation and our previously 
published branching fractions and Af q q is the num- 
ber of background events estimated from sideband data 
(M bc < 5.26 GeV/c 2 ). This TZ requirement retains 74%, 
68%, 74%, and 57% of the signal while rejecting 91%, 
94%, 91%, and 95% of the continuum background for 
the uiK + , u!tt + , cuK°, and lutt modes, respectively. 

Simulation studies indicate small backgrounds from 
generic b — > c transitions in the charged B modes; they 
are found to be negligible for the neutral modes. Two 
additional backgrounds have also been considered: re- 



flections of B — > lott + decays due to 7r — > K misidcnti- 
fication, and feed-down from other charmless B decays, 
predominantly B -> uA"*(892) and B -> ^(770). We 
include these three components in the fit used to extract 
the signal. 

The signal yields and partial rate asymmetries are ob- 
tained using an extended unbinned maximum likelihood 
(ML) fit for M bc and AE. The likelihood is defined as 



(2) 



where i is the identifier of the i-th event, Pj (AA>ci, AEi) is 
the PDF of Mbc and AE, Afj is the number of events for 
category j, which corresponds to either signal, qq back- 
ground, reflections due to K-tt misidentification, or BB 
background (b — > c and charmless). There is no reflec- 
tion component or b — > c background for the neutral B 
modes. 

The signal distribution in Mb c is parameterized by a 
Gaussian function centered near the mass of the B me- 



son, while the Crystal Ball line shape [17j is used to model 
the AE distribution. Both functions are calibrated with 
large control samples of B — > D°ir, D° — > K + ir~iT , 
and B — > rj'K, rj' — > jjtt + tt^, r) —* 77 decays. The 
continuum PDF is the product of a first-order polyno- 
mial for AS and an ARGUS function (T^ for M hc . All 
parameters of the continuum PDF are allowed to float 
in the fit except for the endpoint of the ARGUS func- 
tion, which is determined using the B — > D°ir sample. 
Other background PDFs are modelled by a smoothed 
two-dimensional M^ c — AE function obtained from MC 
simulation. The yields of the continuum and BB back- 
grounds are floated in the fit, except for the reflection 
component. The normalizations of the reflection compo- 
nents are fixed to expectations based on the B + — > u>K + 
and B + — ► uiir + branching fractions, as well as K + <-> 7r + 
fake rates. The reflection yields arc first calculated ac- 
cording to the assumed u>K + /?t + branching fractions and 
then re-estimated based on our measured results from the 
fit. After a few iterations, the reflection yields converge. 
For the charged B modes, we perform independent fits 
to the B + and B~ samples in order to extract Acp ac- 
cording to Eq. 1. 

Figure 1 shows the data and fit results for each mode. 
Table U lists the resulting event yields, efficiencies, cor- 
responding branching fractions and Acp- These new re- 
sults are consistent with those measured by BaBar Q. 
The branching fractions are computed as the sum of the 
yields divided by the corresponding efficiencies in each 
dataset and the total number of B mesons. The recon- 
struction efficiency for each mode includes all interme- 
diate branching fractions and is defined as the fraction 
of the signal yield remaining after all selection criteria, 
where the yield is determined by performing the un- 
binned maximum likelihood fit on a MC sample. Small 



correction factors are included to account for differences 
in the reconstruction efficiency between MC and data. 
The correction factor is obtained as a result of the relative 
7T reconstruction (1.00 for Set I and 0.96 for Set II) using 
the inclusive D° — » K + ir~ and D° — > K + ir~ir° samples. 
The corresponding systematic error is 4%. We determine 
the identification efficiency of charged kaons and pions by 
studying an inclusive D*~ — > D°ir~ , D° — > K + tt~ sam- 
ple. This selection criterion leads to a correction factor 
within the range of 0.90-0.96 and a corresponding sys- 
tematic error of 1.7% for charged modes and 0.9% for 
neutral modes. 



We define the significance as yj — 21n(£o / Anax) , where 
£ ma x is the maximum likelihood from the fit when the 
signal yields are floated, and Co is obtained when the 
signal yields are set to zero with the remaining fit pa- 
rameters still floating. To verify that signal candidates 
indeed arise from B — > cuK or B — > tun, we in turn change 
the criterion on M^^+^-^o or | cos 0£ el \, and repeatedly 
perform ML fits for M^ c and AE to obtain yields in each 
bin of A^^tt+tt-ttO or | cos0£ ol | shown in Figure 2, where 
&£ el denotes the angle between the B flight direction and 
the normal to the decay plane in the to rest frame. The 
90% C.L. upper limit for B(B° — ► ojit ) is calculated 
using the procedure of Ref. 19(, which is based on the 
Fcldman-Cousins method [20|. This procedure accounts 
for systematic uncertainties (discussed below). The es- 
timated background (observed count) is 57 (65) events 
inside the M^ c — AE signal box. 

The main systematic uncertainties for the branching 
fraction measurements are listed in Table [TTJ These are 
evaluated as follows: charged tracking efficiency (1.0% 
per track) from partially reconstructed D* + — > D°tt + , 
D° — > K~it + ; Kg reconstruction (4.9%) from D + — > 
Kgir + , D + — > A'~7r + 7r + ; uj mass resolution (3.0%); the 
requirement on 1Z (3.0%) from B~ — > D q tt~, D° — > 
K-TT+; MC statistics (1.1-2.4)%, and the number of BB 
events in the data samples (1.3%). For the Acp measure- 
ment, the dominant systematic uncertainty comes from 
the fit parameters. Systematic uncertainties due to fit 
parameters are evaluated by varying each parameter by 
±1(7 and adding in quadrature the resulting deviations 
from the central value. The total systematic uncertainty 
is the quadratic sum of all the above contributions. When 
varying each fit parameter by Ha one at a time, the sig- 
nificance is recalculated and the lowest value is chosen 
as the final significance (shown in Table [l| including the 
systematic uncertainty. 

While the branching fractions for B + — > uiK + , B° — > 
ujK°, and B + — > lott + are of comparable size, that for 
B° — » ujtt is much smaller as compared to that for B + — > 
cj7r + . This is in agreement with both QCDF and PQCD 
pictures, which predict that only color-suppressed tree 



and penguin diagrams contribute to B 
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0,0, |. 



TABLE I: Signal efficiency (e), signal yield (j\f s ), significance including systematic error (E s i g .), branching fraction (B), the 
90% C.L. upper limit (U.L.) and partial rate asymmetry (Acp) for the B + — > coh + decays. The first and second errors are 
statistical and systematic, respectively. 



Mode 


e (%) Ms 


S s ig. 


B (W- 6 ) 


U.L. (HT 6 ) 


Acp 


ujK+ 


8.27 ±0.17 259.5iig;t 


19.5 


8.1 ±0.6 ±0.6 




0.05±g;g| ± 0.01 




8.43 ± 0.17 224.8±i9.3 


17.1 


6.9 ±0.6 ±0.5 




-0.02 ±0.09 ±0.01 


loK° 


2.50 ±0.09 41.5l?;g 


9.3 


4.4l°;* ±0.4 






UJTt° 


3.80 ±0.12 5.9t£? 


1.5 


o.5i°:*±o.i 


2.0 





150 - 




TABLE II: systematic uncertainties for ujK and ujtt (%). 



> 5.2 5.225 5.25 5.275 5.3 

m M bc (GeV/c 2 ) 



0.2 

AE (GeV) 



FIG. 1: Projections of fit results on Mbc (within the region of 
\AE\ < 0.1 GeV, -0.15 GeV < AE < 0.1 GeV for ujtt ) and 
AE (within the region of 5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 ) 
for uiK + , un + , luK°, and loty". Open histograms are data, 
solid curves are the fit functions, solid-dotted lines represent 
the signals, dashed lines show the sum of qq continuum and 
b — > c, dotted lines are due to other charmless B decays, and 
cross-hatched areas represent reflections due to K-ir misiden- 
tification. 



Mode 


loK + 


LJTT + 


ujK 




UJTT 


track reconstruction 


3.0 


3.0 


4.0 


2.0 


1Z requirement 


2.5 


2.6 


2.5 


3.3 


particle identification 


1.7 


1.7 


0.9 


0.9 


7T° reconstruction 


4.0 


4.0 


4.0 


8.0 


Kg reconstruction 






4.9 




uj mass window 


3.0 


3.0 


3.0 


3.0 


MC statistics 


1.1 


1.2 


1.8 


2.4 


signal PDF 


+0.1 
-0.0 


< 0.1 


±0.5 


+ 1.2 
-0.6 


uiK + /um + feed-across background 


+0.3 
-0.1 


±0.3 






N B B 


1.3 


1.3 


1.3 


1.3 


Total 


6.8 


6.8 


8.8 


9.9 



We also determine the partial width ratio 



1.7 ± 



0.3(stat.) ±0.1(sys.), (3) 



where the ratio of B meson lifetimes 
be 1.076 ± 0.008 



12] 



•b + 



is taken to 



The systematic error of the par- 
tial width ratio is reduced because of the cancellation 
of several common systematic errors. No evidence for 
direct CP violation is found for either B + — > loK + or 
B + — » lott + . Our results on branching fractions and 
CP asymmetries are consistent with the recent measure- 
ments Q from the BaBar collaboration. 
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FIG. 2: Background-subtracted one-dimentional projections 
are obtained from data fit results in bins of fit yields in 
(left) the invariant mass M^^^+^-^o and (right) the ab- 
solute value of the cosine of the helicity angle | cos 9^ A | for 
luK + , uj-k + , and ujK° candidate events. Points are data, and 
histograms are from MC expectations for B — > uih. 
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